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Summary. For a given orthonormal basis (/„) on a probability measure space, 
we want to describe all Markov operators which have the /„ as eigenvectors. We 
introduce for that what we call the hypergroup property. We study this property in 
three different cases. 

On finite sets, this property appears as the dual of the GKS property linked 
with correlation inequalities in statistical mechanics. The representation theory of 
groups provides generic examples where these two properties are satisfied, although 
this group structure is not necessary in general. 

The hypergroup property also holds for Sturm-Liouville bases associated with log- 
concave symmetric measures on a compact interval, as stated in Achour-Trimeche's 
theorem. We give some criteria to relax this symmetry condition in view of extensions 
to a more general context. 

In the case of Jacobi polynomials with non-symmetric parameters, the hyper- 
group property is nothing else than Gasper's theorem. The proof we present is based 
on a natural interpretation of these polynomials as harmonic functions and is related 
to analysis on spheres. The proof relies on the representation of the polynomials as 
the moments of a complex variable. 
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1 Introduction 

In a number of situations, Markov operators appear to be a wonderful tool 
to provide useful information on a given measured space. Let us for example 

mention heat kernel methods to prove functional inequalities like Sobolcv or 
Log-Sobolev inequalities, or Cauchy kernels to prove boundness results on 
Riesz transforms in LP. Heat kernels are widely used in Riemannian geometry 
and statistical mechanics, while Poisson, Cauchy and other kernels had been 
proved useful in other contexts related to classical harmonic analysis (see 
[9, 5, 6, 7, 11, 10, 20, 28, 29, 2, 38] for example to see the action of different 
families of semigroups in various contexts). 

It seems therefore interesting to describe all Markov kernels associated 
with a given structure. In what follows, we shall consider a probability space 
{E, E, /i) on which is given an orthonormal L'^^fJ.) basis J- = (/o, . . . , /«, . . .), 
where we impose /o = 1. Such a basis shall be called a unitary orthonormal 
basis (UOB in short). In general our basis shall be real, but we do not 
exclude to consider complex bases. 

We may then try to describe the Markov operators defined from a family 
of probability measures k{x,dy) by 

f^K{f)= [ f{y) k{x,dy) 

JE 

which are symmetric in L'^ip) and have the functions /„ as eigenvectors. In 
other words, we want to define the linear operator K from 
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and try to describe for which sequences (A„) this operator is a Markov kernel. 
We shall call these sequences Markov sequences (MS's in short) associated 
with the UOB T. 

For a general basis JF, this is quite impossible. But many bases which 
appear in natural examples have a special property, which we call the hy- 
pergroup property, under which one is able to describe all Markov sequences 
associated with the UOB J^. 

This expository paper is not intended to be a complete account of the 
general theory of hypergroups, for which we may for example refer to the 
complete treatise [13]. In fact, we just extracted from this theory what is 
useful for our purpose. More precisely, we concentrated on the fundamental 
aspect which wc arc interested in, that is the possibility of describing all 
Markov kernels associated with our basis J^. 

The paper is organized as follows. 

In the first part, we present the case of finite sets, where the hypergroup 
property appears as the dual property of a more natural condition on the 
basis namely the positivity of the multiplication coefficients. This prop- 
erty is called the GKS property in [8] because of its links with some famous 
correlation inequalities in statistical mechanics, and we keep this notation. 
These correlation inequalities did in fact motivate our efforts in this direction 
(see Paragraph 2.6). Many examples come from the representation theory of 
groups, but we propose a systematic exposition. The hypergroup property 
provides a convolution operation on the set of probability measures on the 
space, and the Markov kernels may be represented as the convolution with 
some given measure. In this situation, the hypergroup property appears as a 
special property of an orthogonal matrix. We shall see that in fact there arc 
many situations where no group structure holds and where nevertheless the 
hypergroup property holds. 

The second part is devoted to the presentation of Achour-Trimeche's the- 
orem, which states this property for the basis of eigenvectors of a Sturm- 
Liouvillc operator with Numen boundary conditions, associated with a log- 
concave symmetric measure on a compact interval. The original Achour- 
Trimcchc's theorem was not stated exactly in the same way (sec [1]), but 
his argument carries over very easily to our context. We give a complete proof 
of this result, since to our knowledge this proof was never published. We 
tried to relax the symmetry condition on the measure, and provided various 
rather technical extensions of the theorem. But wc did not succeed to extend 
Achour-Trimeche's theorem to a wider class of measures which would include 
the case of Gasper's theorem on Jacobi polynomials studied in the next chap- 
ter. The real motivation of this effort is that in general, log-concave measures 
on R or on an interval serve as a baby model in Riemannian geometry for 
manifolds with non-negative Ricci curvature. Unfortunately, the symmetry 
condition does not seem to have any natural interpretation. 
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In the last part, we present Gasper's theorem, which states hypergroup 
property for the Jacobi polynomials. We present a proof which relics on ge- 
ometric considerations on the spheres when the parameters are integers, and 
which easily extends to the general case. We follow Koornwinder's proof of the 
result, and give a natural interpretation of Koornwinder's formula (Lemma 
5.2) which represents those polynomials as the moments of some complex ran- 
dom variable. We found after the redaction of this part that the interpretation 
of Jacobi polynomials as harmonic functions was already known from special- 
ists (see [14, 23, 35]) but it seems that it was not directly used to prove this 
integral representation formula. We hope that this simple interpretation may 
provide other examples for similar representation in other contexts. 

2 The finite case 

In this section, we restrict ourselves to the case of finite sets, since in this 
context most of the ideas underlying the general setting are present, and we 
so avoid the analytic complexity of the more general cases that we shall study 
later on. 

2.1 The GKS property 

In what follows, we assume that our space is a finite set 

E = {xq, . . .,Xn}, 

endowed with a probability measure 

fj,= (/i(xo),...,M(a;„)). 

We denote by L'^{fJ,) the space of real functions on E, and we assume the 
existence of a real basis 

^ = (/O) • • • ) fn) 

with /o = 1. We suppose here that for any x G E, fi{x) > 0. 
We shall write (/) for / / and (/, g) for / fg dji. 

The algebra structure of the set of functions is reflected in the multiplica- 
tion tensor (a^fe) for which 

fifi — ^ ] dijkfk- 
k 

We therefore have 

aijk = ififjfk), 

and we see that the tensor (aijk) is symmetric in (ijk). It has also another 
property which reflects the fact that the multiplication is associative. 
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Definition 2.1. We shall say that T has the property GKS if all coefficients 
(aijk) are non-negative. 

This notation comes from the GKS inequahty in statistical mechanics that 
we shall describe at the end of this section. 
Observe that = • 

Many natural bases !F share this property. For example, consider the hy- 
percube E = {—1,1}^, with the uniform measure on it. Let denote the 
i-th coordinates 

LJ = (Wl, . . .,U!n) 1-^ U>i, 

and, for Ac {1,...,N} 

UJA = J]^ Wi, UliD = 1. 

Then, 

T={u;^uJA,Ac{l,...,N}} 
is a UOB of {E,ii). Since 

it has the GKS property. 

(Wc shall see later that this is a special case of a generic situation in finite 
groups). 

Although we are here mainly interested in the case of a real basis, there are 
many natural complex GKS bases, issued in general from the representation 
theory of finite groups (see Paragraph 2.5 later). If the basis is complex, we 
shall still require that the multiplication coefficients are non-negative real 
numbers, which means that 

O-ijk = ififjfk) > 0, 

for any {i,j,k). 

In what follows, we only consider real GKS bases, although the next result 
remains probably true in the complex setting. 

Proposition 2.1. If a UOB has the GKS property, then there exists a unique 
point xq on which every fi{xo) is maximal. Moreover, for any i and any x, 
\fi\ {x) < fi{xo), and at this point xo, ij,{xo) is minimal. 

Proof. Let us say that a function / : _B M is GKS if for any i = 0, . . . , n, 
iffi) > 0. In other words, / is written with non-negative coefficients in the 
basis !F. 

We shall say that a set X C -E is GKS if Ix is a GKS function. We shall 
say that a point x is GKS if {x} is a GKS set. 
We shall see that there is only one GKS point. 

Remark first that the STim of two GKS functions is GKS and that, thanks 
to the GKS property of J^, the product of two GKS functions is GKS. Moreover 
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a limit of GKS functions is GKS. Observe also that a GKS function has always 
a non-negative integral with respect to fi since /o = 1 • 

Let us consider a non-zero GKS function / and consider m = maXx^E |/|- 
We see first that {/ = m} 7^ 0. 

For this, assume the contrary, that is that / = — m on |/| = m. Since m > 
0, we see that /^^+^/m^^"'"^ is a GKS function, and converges to — 
Since a GKS function has a non-negative integral, this is impossible. Using 
the same argument, we see that 



is a GKS function. 

Therefore, the set {/ = m} is a GKS set, and there are non-trivial GKS 
sets. 

Moreover, for any GKS function, the set {/ = max |/|} is GKS. 

Let El be a nonempty GKS set, minimal for the inclusion. Then, for any 
GKS function f , g = Ie^J is GKS. If g is not 0, then its maximum is attained 
on a subset E2 of Ei which is again GKS. Since Ei is minimal, we have 
E2 = El. 

Therefore, for any GKS function, its restriction g to Ei is either on Ei 
or constant (and equal to the maximum of 9). In any case, / is constant on 
El. Since this applies to every function in J^, and since is a basis, every 
function is constant on Ei and therefore Ei is reduced to a single point {xq}- 

The same proof shows that any GKS set contains a GKS point. 

For a GKS point Xq, fi{xQ) = {fila:o) / h{xq) > 0. 

Then, consider two distinct points xq and Xi , and write 



Suppose then that xq and xi are GKS points. In the previous sum, all 
coefficients are non-negative. Therefore, for any k) 





Writing the product, we see that 




with the multiplication coefficients aijk. 

So we see that for any pair of distinct points, and for any k, 



^fi{xo)fj{xi)aijk = 0. 



aijkfi{xo)fj{xi) = 0. 
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If we apply that with i = and j = k, we see that fj{xi) = 0, for any i. This 
is impossible since then any function would take the value in Xi. So there 
is a unique GKS point. 

Let xo be this unique GKS point. Any GKS set contains xq- Since for any 
GKS function, the set where / is maximum is GKS, any GKS function attains 
its maximum at xq. 

It remains to show that /x is minimal at xq. For this, observe that for any 
point X, the function 

fx = ^J-ix)lxa - K^o)'ix 

is GKS (this comes from the fact that each /j is maximal at Xq). Therefore, 
the maximum value of is attained in xq, which gives the result. □ 

2.2 Orthogonal matrix representation 

Consider the matrix 

{Oij) = (^y/jMx^fj{xi)j , 

we see easily that the matrix (Oij) is a (n+ 1) x (n+ 1) orthogonal matrix with 
positive first column. Conversely, any such matrix may be associated with a 
UOB on a finite set with measure /x given by 

nix,) = OIq. 

Therefore, there is a one to one correspondence between the set of orthogonal 
matrices with positive first column, and the set of finite probability spaces, 
whose probability has everywhere positive weight, endowed with a UOB. (In 
fact, this is not completely true, since we would not distinguish between bases 
given in different orders, provided that the first element is 1, which identifies 
the set of UOBs with a quotient of a the set of orthogonal matrices through 
a permutation of rows and columns.) 

The GKS property may be translated into the following property on such 
an orthogonal matrix: 

The transposed of an orthogonal matrix is orthogonal, and we just saw that 
an orthogonal matrix which has the GKS property also has a non-negative 
row (corresponding to the row where the first column is minimal according 
to Proposition 2.1). We may of course rearrange the labelling of the points 
in such a way that this row is the first one. Then, the situation is completely 
symmetric. 

We shall then consider the squares of terms in the first row as a probability 
measure on the dual set {0, 1, . . . , n}: 



v{i) = i^{xo)f^{xo). 
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Thanks to the fact that the functions /j are maximal at xq and that this 
maximum must be larger than 1 (since J ff diJ,= 1), the dual measure is also 
minimum at 0. 

As an application, we have the following. 

Proposition 2.2. If a real GKS basis exists for the uniform measure on some 
finite set E , then the cardinal of E must he 2^ for some k, and this basis is 
the canonical basis {u)a) of the characters of the group (Z/2Z)'^. 

To sec this, we first observe that the dual measure is uniform too. In fact, 
for the matrix (Oy), we have Ooo = where n + 1 is the number of 

points in the space, and since the first row is positive and has minimum value 
1 / -v/rT+T, it must be constant since the sum of the squares of its coefficients 
is 1. 

Now, if we multiply the matrix O by y/rT+T, then we see that in each 
column, the maximum value of the coefficients is attained on the first row 
and is equal to 1. Since the sum of all the squares of the coefficients in a given 
column must add to n + 1, this shows that in any column, the coefficients 
must take only the values ±1. 

Those matrices (with entries ±1 and orthogonal lines) are called Hadamard 
matrices (cf [30, 37, 33]). If n + 1 = 2*^ for some k, such matrices are given 
by the basis (lua) on { — 1, 1}*^ and satisfy the GKS property. It is known 
that the order of a Hadamard matrix must be 1, 2, or 4fc, and it is an open 
problem to find such matrices for all k (the lowest k for which no Hadamard 
matrix of order 4k is known is fc = 167 since [33]). Nevertheless, the following 
proposition will prove the result of Proposition 2.2. 

Proposition 2.3. // a Hadamard matrix has the GKS property, then it must 
be of order 2^ for some k, and, up to permutation, it is the matrix of the 
canonical basis {coa) of the group {—1,1}'^. 

Proof. The case of a set of size 2 is trivial, and we therefore assume that the 
size of the matrix is at least 3. To fix the idea, consider a matrix M of order 
n + 1 with entries ±1, with orthogonal columns. Call fi the column vectors 
and fi = fi/Vn+ 1 the normahzed ones, so that M = (.^(■^j))o<i j<„- 
usual, let us denote by (/) the mean value of a function / with respect to the 
normalized uniform measure. We suppose, which is possible up to reordering, 
that all the entries of the first line and the first column are +1. The GKS 
property says that 

{f^fJfk) > 0, 

for any k). 

First, since {fi) = for any i > 1, there must be as many 1 and —1 in each 
column, and therefore n + 1 is even. Let A be the set of points where /i = 1. 
Let i > 2, and let q be the number of points in x G A such that fi{x) = 1, 
and r be the number of points in A" where fi{x) = 1. Writing (/j) = and 
(/i/i) = 0, we get 
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q + r=p, q = r, 
which shows that p is even and also that, when i>2, 

{tAf^) = 0. 

This is the generic argument which shows that Hadamard matrices have order 
4fc. We shall now make use of the GKS property. Write 

We have ^ 

an = {^Afifi) = — —ri'^Afifi), 
n + 1 

and therefore the matrix A — {an) is symmetric. As it is the matrix of a 
projector, its eigenvalue are or 1. The dimension of tlic; cigcnspacc associated 
to 1 is p, since the eigenspace is generated by (l^^,}, a; G A). Also, A is a GKS 
set because it is the set where /i attains it's maximum (cf proof of Proposition 
2.1, Page 6). This implies that all the entries of A are not negative. Moreover, 
if we look at the values of the functions at xq, it holds Ylij ^ij = 1. 

Therefore, the matrix A is Markovian. There are no transitory points since 
A is symmetric. The number of recurrence classes for such a matrix is the 
multiplicity of 1 as eigenvector, herep, so there are exactly p recurrence classes, 
and no recurrence class is reduced to a single point, since an = 1/2 < 1. 
So every recurrence class has exactly two points. For example, {0, 1} form 
a recurrence class. On each line of the matrix, there are exactly two places 
where a^j 0. The values of those entries are then 1/2, since an = 1/2. If 
we choose two distinct indices i and j in two different recurrent classes, then 
aij = 0. This means that {lAfifj) =0. 

Choose now one index in every recurrence class (say the even indices to 
fix the ideas, which is possible up to reordering of the columns). Then, the 
functions gi = lAf2i form an orthogonal Hadamard GKS matrix of order 
p = (n + l)/2, and we may now use induction to see that the order must be 
2*^ for some k. 

To see that the unique basis such basis in dimension 2'' is given by the 
canonical basis, it is enough to observe that if i and a{i) are in the same 
recurrence class for the matrix A, then we have 

Ia/j = -^{fi+ fcr(i)) = lA/<T(i), 

and also 

Then an easy induction leads to the result. □ 
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Nevertheless, unlike real bases, we shall see in Paragraph 2.5 that there 
always exists a complex GKS basis on any finite set with uniform measure. 

On two points, an easy computation shows that a two dimensional orthog- 
onal matrix having the GKS property must be 



with 6 e [7r/4, 7r/2). In fact, in dimension 2, given the measure (giving two 
distinct positive masses on the two points), there arc exactly two UOB. and 
only one such that the unique non-constant function is maximal on the point 
with minimal mass (a necessary condition to have the GKS property as we 
saw). In this situation, any GKS matrix is symmetric, and the set of orthogonal 
matrices having the GKS property is connected. 

On three points, the situation is more complicated. We saw for example 
that there arc no real GKS basis when the measure is uniform. The set of real 
GKS UOBs on three points is connected, and one may see that the maximum 
value of n{xo) for which there exists a real GKS basis is ij,{xq) = 1/4, and in 
this case the probability measure is (1/4, 1/4, 1/2) and the unique GKS basis 
is obtained taking the real parts of the characters in the group Z/4Z (we 
shall see later in Paragraph 2.5 how to associate complex or real GKS bases 
with any finite group). There is also a complex GKS basis with the uniform 
measure (the characters of the group Z/3Z). 

2.3 The hypergroup property 

As we saw before, when we have the GKS property, the situation is completely 
symmetric and we may consider the dual property. This is the hypergroup 
property. 

Definition 2.2. We shall say that the UOB T of L^{^) satisfies the hyper- 
group property (HGP in short) at point Xq if for any {x,y,z) G , 



Of course, this supposes that for any i, fi{xo) ^ 0. 

We do not require the GKS property to hold in the definition of the HGP 
property. We shall see later (at the end of Paragraph 2.5, Page 18) that we 

can find bases with HGP property and not GKS property, and the reverse. 
On the other hand, we shall see that in groups, the natural basis always share 
both property. 

Remark. Since we may always change for any i > 1 fi G into — /j, and still 

get a UOB, and that this operation does not change the hypergroup property, 
we see that we may always assume that fi{xo) > 0. 



cos 6 sin 
sin 6 — cos 9 
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To see the duality with the previous situation, let us enumerate the 
points in E starting from Xq, and recall our orthogonal matrix O with 

Then, under the GKS condition, O has non-negative first line and first 
column. Recall that the GKS property may be written as 



Notice also that if both properties occur, then the point xo must be the 
unique point where all functions fi are non- negative, and where all and |/,| 
are maximal. 

Prom the symmetry of the situation, we may consider the functions 



to be a UOB on the set {0, . . . , n} endowed with the measure = O^^. 

From this we deduce that a basis has the HGP property if and only if this 
new basis has the GKS property, where of course the point plays the role of 
the point Xq in the previous paragraph. Therefore, the gx(i) are maximal at 
i = 0, and moreover Ojo is minimal at i = 0. This means that, for the HGP 
property also, one has 

Proposition 2.4. // the UOB T has the HGP property at point xq, then 

Vi,x, \fi{x)\ < fi{xo); yx, ii{x) > n{xo). 

(Recall that we assume that fi{xo) > 0.) 

We may reformulate the HGP property in the following way, which we 
shall use later in a different context, since it is more tractable. 

Proposition 2.5. The UOB T has the HGP property if and only if there is 

a probability kernel k{x, y, dz) such that, for any i = 0, . . . ,n 



Proof. The proof is straightforward. If the hypergroup property holds, then 
the kernel 




while the HGP property writes 
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is a probability kernel satisfying our conditions. 

On the other hand, if such a probability kernel k{x,y,dz) exists, then 
writing 

k{x, y, dz) = Ki{x, y, z)n{dz) 



and 



one sees that 



Ki{x,y,z) = ^ai{x,y)fi{z), 

i 

Mx)My) 



ai{x,y) 



fi{Xo) 



□ 

The link with the Markov operators is the following. 

Definition 2.3. A Markov operator is just an operator K which satisfies 
K{1) = 1 and which preserves non-negative functions. 

Given a sequence X = (Aj, i = 0, . . . ,n), we define the associated linear 
operator Kx by 

Kxifi) = \h. 

A Markov sequence (MS in short) is a sequence A = (A,) such that the asso- 
ciated operator K\ is a Markov operator. 

Remark that for any MS A, one has Aq = 1. Remark also that 

Kx{f){x) = J f{y) kx{x,dy), 

where 

k\{x,dy) = i^^\ifi{x)fi{y)^ ii.{dy). 

Therefore, the set of MS's is just the set of sequences A such that the matrices 
k\{x,y) axe Markov matrices. 

The HGP property asserts that, for any x, the sequence 



X{x) 



M^O) J i=0,...,n 



is a MS. 

It is quite standard to see that any eigenvalue Aj of a Markov operator 
must satisfy |Ai| < 1. The set of Markov sequences is a convex compact set, 
which is stable luidcr pointwisc niidtiplication. Now the main interest of this 
property relies in the following theorem. 
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Theorem 2.1. // the basis T has the HGP property, then the sequences 



Xi{x) = 



fi{xo) 



are the extremal points in the convex set of all Markov sequences. 
More precisely, every Markov sequence may he written uniquely as 

. f Mx) , , . 

for some probability measure v on E. Conversely, every probability measure 
can be associated in the same way with a Markov sequence. 

The main interest of this result is that there exist numerous natural 
bases with the HGP property, as we shall see later. 

Proof. The representation formula is straightforward. Indeed, writing ^{dy) = 
kx{xo, dy), one has 

Kfi{xo) = j fi{y) v{dy), 

which gives the representation. 

Prom this, it is easy to see that if the sequences {\i{x))i are Markov 
sequences, then they are extremal. Indeed any representation 

with MS's and leads to 

fi{x) = j fi{y) v{dy), 

with 

v{dy) ^ eX^i {xo,dy) + (1 - 9)K),. {xq, dy). 
From this we deduce that for any function / 



fix) = I f{y) u{dy), 



and therefore ly ~ S^, which gives the extremality. □ 

Remark. Remark that the representation formula is still true for any MS 
when the basis does not verify the HGP property. As we can always embed 
the convex set of Markov sequences in the n-dimcnsional affinc space Hi ~ 
{(Aj), Ao = 1}, the latter fact means that this set is actually contained in the 
n-simplex generated by the n + 1 points A^/^° = {fi{x)/fi{xo)).. 
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Then, when the hypergroup property holds, the set of Markov sequences is 
a n-simplcx and the representation of a point in this set as afRne combination 
of extremal points is unique. 

We may ask for which kind of L"^ basis on a finite space this still happens. 
It is quite clear that the cardinal of the set of extremal Markov sequences is 
finite. Indeed, the set of Markov sequences is delimited by a finite number 
of (n — l)-hyperplanes in Hi. Namely, for any pair {x,y) of points in E, one 
considers the half space defined by {(A, F^'^) > 0}, where F^'^ = {fi{x)fi{y)). 
Then the set of Markov sequences is the intersection of all these half spaces. 
Therefore, every extremal point lies in the finite set Ei of possible intersections 
of n hyperplancs H^^y = Hi D {(A,F^'^) = 0}. Now, consider any point xo 
such that for any index i, fiixo) ^ 0. The point \^/^° = {fi{x)/fi{xQ)). 
belongs to iJ^"'" for any y ^ x, thanks to the orthogonality relations of the 
basis. Therefore, those points \^/^« belong to the set Ei. When xo is fixed 
and X varies in E, those points describe a simplex Sx^ for which we know that 
every Markov sequence belongs to it. The hypergroup property holds at some 
point xo exactly when no other point in Ei lie in the interior of Sx„ ■ 

On three points, one may check directly that the hypergroup property 
holds at some point xq exactly when the set of Markov sequences is a simplex 
(that means no other simplex is possible than the simpliccs Sx, x E E). We 
may wonder if this situation is general, that is if the hypergroup property is 
equivalent to the fact that the set of Markov sequences is a simplex. 



2.4 Markov operators as convolutions 

When the hypergroup property holds, we may introduce a convolution on the 
space of measures. 

Indeed, consider the kernel 

We observe that, for any (x, y) 

j k{x,y,z) ii{dz) = 1, 

and therefore the measures 

fj,x,y{dz) = k{x,y,z) iJ,{dz) 

are probability measures. 

We may decide that the convolution is defined from this kernel by 

Sx * Sy ^ Hx.y, 



and extending it to any measure by bilinearity. 
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Moreover, we extend the convolution to functions by identifying a function 
/ with the measure fd^. This gives 

/ * 9{z) = j f{x)g{y)k{x, y, z) dii{x)dii{y). 

Observe that 

and that this property again completely determines the convolution. 
It is easy to verify that this convolution is commutative and that 5xo*v = v 
for any u. Moreover, if an operator K satisfies K{fi) = Xifi, then 

K{f*g)=K{f)*g = f*Kig), 

as may be verified directly when / = and g = fj using (2). 

On the other hand, if is a probability measure, then the operator K^{f) = 
f * u is a Markov operator which satisfies 

KM^) = hh, 

with 

^ ^ //i dv 
fi{xo) ' 

This is straightforward using (2) if we write v{dx) = h{x)fj.{dx) and the de- 
composition of h along the basis J^. 

Therefore, if -fT is a Markov operator, then we have 

This representation is exactly the representation of Markov sequences, with 
u = K{Sxg), and every Markov operator Kx may be defined from K\{f) = 
f * v, for some probability measure v. 



2.5 The case of finite groups 

Many natural examples of finite sets endowed with a probability measure and 
a UOB which satisfies both GKS and HGP properties come from finite groups. 

Since perhaps not every rc^adca' of these notes is familiar with this setting, 
let us summarize briefly the basic elements of the analysis on groups. We refer 
to [22] or [31] for more details. 

Given a finite group G, one may consider linear representations p : G ^ 
C/(V), that is group homomorphisms between G and some UiV), for some 
finite dimensional Hermitian space V (where U{V) denotes the unitary group 
of V). Such a representation is irreducible if there is no non-trivial proper 
subspace of V which is invariant under p{G). Any representation may be split 
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into a sum of irreducible representations, acting on orthogonal subspaces of V. 
Two representations {pi, Vi) and {p2, V2) arc equivalent if there exists a linear 
unitary isomorphism h : Vi 1-^ V2 such that pi{g) = h~^p2{g)h for any g £ G. 
There are only a finite number of non-equivalent irreducible representations, 
that we denote {pi, Vi), i £ I ~ {0, . . . , n}. 

Let G the set of the equivalence classes of G under the conjugacy relation 
{gi is conjugate to 52 means 52 = g~^gig for some g G G). We endow G with 
the probability v which is the image measure of the uniform measure on G, 
which means that the measure of any class is proportional to the number of 
points in this class. A function on G which is constant on conjugacy classes 
(that we call a class function) can be seen as a function on G. It is just a 
function which is stable under conjugacy. 

For any irreducible representation {pi,Vi), let us define the function Xi on 
G by Xi{g) = trace This is a class function, that is to say constant 

on any conjugacy class. The function Xi is called the character of the repre- 
sentation. By convention, we take xo = 1; that is the trace of the constant 
representation into the space V = {C}. 

Proposition 2.6. The set {xi- i £ 1} is a (complex) UOB for {G,u). More- 
over, it has the GKS and HGP properties. 

Proof. We shall not enter in the details here. We refer to any introduction book 
on the representation theory of finite groups for the first fact. We shall detail 
a bit more the HGP and GKS properties, which are perhaps less standard. 

For the GKS property, for any pair of irreducible representations {pi,Vi) 
and {pj, Vj), one may consider the representation pi® pj in the tensor product 
Vi (g) . If we split this representation into irreducible representations and take 
the trace, and if we notice that trace [pi P2) = trace (pi)trace (^2), then we 
get that 

XiXj ~ ^ ] ^ijfcXfc; 

k 

where rriijk is the number of times that the representation p}. appears in this 
decomposition. Here we may see that not only the multiplication coefficients 
are non-negative, but they are integers. 

We shall see next that this basis has the HGP property at the point xq = e 
(which forms a conjugacy class by itself). For that, we require a bit more 
material. 

First define the convolution on the group G itself by 

I I g'GG 

The Fourier transform is defined on the set I of irreducible representation 

as 



The Hypergroup Property and Representation of Markov Kernels 17 

(It takes values in the set of linear operators on V^.) 
One has an inversion formula 



Hg) = j^^Y^ditYcice {pi{g 



where di is the dimension of Vi (the degree of the representation) . 
One has 

(^* v]r= ^V', 

and 

Now, the convolution of two class functions is again a class function, as 
seen directly from the definition. 

We want to show that this convolution is exactly the convolution that we 
defined in the previous section from the HGP property, that is 

For that, wc look at the Fourier transform and the result is straightforward, 
since Xi(e) = d^. 

This convolution is then the convolution defined from the Xi, and we have 

5x*5y = k{x, y, z)dii{z), 

where 

If \ \^ Xi{x)Xi{y)Xi{z) 
k{x,y,z) = j: . 

Since by construction in this case the convolution of two probability mea- 
sures is a probability measure, the kernel k{x,y,z) is non-negative, which 
proves the HGP property. 

Observe that here the kernel k{x, y, z) has a simple interpretation. Given 
3 classes {x,y,z), then 

k{x,y,z) = -^^m{x,y,z), 

m \y\ 

where m{x, y, z) is, for any point g ^ z, the number of ways of writing g = gig2 
with gi x and g2 S y, this number being independent of the choice of g z. 
□ 

If wc want to stick to real bases as we did before (and as shall do in the 
next chapters), we may restrict ourselves to real groups (that is groups where 
g and g~^ are always in the same class), or we may agglomerate the class of 
g with the class of .9^^. Wc get a new probability space, where the functions 
3?(Xi) form a UOB which again satisfies the GKS and HGP properties. 
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It is certainly worth noticing that, unhke the convolution on G itself, the 

convolution on G is always commutative. 

Observe that taking the group "L/nL, one gets a complex GKS and HGP 
UOB on the set of finite points with the uniform measure (with fi{x) = 
exp{2iTrlx)), and that the unique real case where the measure is uniform and 
is GKS (the hypercube) is nothing else that the group (Z/2Z)"'. 

Remark. Unlike what happens for finite groups, it is not true in general that 
a basis ^ which has the GKS property has the dual property HGP. This is 
the case on two points spaces, since any orthogonal GKS matrix is symmetric 
(cf Page 10). If we look at the sets with three points, one may construct 
examples of an orthogonal matrix having the GKS property without the HGP 
property (and conversely, of course). In fact, consider an orthogonal matrix 
(Ojj), < i,j < 2, with positive first row and columns. Ooq,Ooi and Oio 
determine entirely the first rows and columns, and then it is easy to see that 
there are only 2 orthogonal matrices with given Ooo,Ooi,Oio- Then, it is 
not hard (using a computer algebra program) to produce orthogonal matrices 
which have the GKS and not the HGP property, or which have neither, or 
both. 



2.6 On the GKS inequalities 

We conclude this section with some remarks on the correlation inequalities in 
statistical mechanics. 

In this context, one is interested in the space of configurations of some 
system. We have a set of positions i G K, K being a finite set, and at each 
point i £ K there is some random variable x, with values in E, where E is 
some finite set, endowed with a probability measure /x. One is then interested 
in the set of configurations, which is equipped with a measure /Ujj, where 

HHidx) = exp(i?(x))— , 

where /xq is the product measure /x'^^ on E^ , H is some function on E^ (the 
Hamiltonian) , and Zh is the normalizing constant. 

One of the basic example of spin systems is when E = {—1,1}, and H = 
cawa, where the functions wa are the canonical GKS basis on {—1, 1}^ 
described before. 

To study such systems (and more precisely their asymptotics when K 
enlarges), one uses some structural inequalities. Wc present here two funda- 
mental such inequalities, known as GKS inequalities, from Griffiths [27], Kelly 
and Sherman [32]. The GKS property for a basis has been introduced in [8], 
in an attempt to generalize the GKS inequality to a more general context. 

The classical GKS inequalities are settled in the context of (Z/2Z)^. As 
before, we say that F is a GKS function if F = Xl^icif /^'^a, where V.A, 

>0. 

Then we have 
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Proposition 2.7. 1. (GKSl inequality). Assume that F and H are GKS. 



2. ( GKS 2 inequality). Assume that F, G and H are GKS functions. Then 



The main advantage of the GKS and HGP properties is that they are 
stable under tensorization. That is, if one considers two sets (i?j,/Ltj) with 
UOB bases J^i {i = 1, 2), then, on the set {Ei x E2, ^1 (Xi /X2) one has a natural 
UOB basis Ti®T2 = {.U ® fj). Then, if both are GKS or HGP, the same 
is true for J^i ® ^2- This is straightforward from the definitions. 

This allows us to consider a set {E, fi) with a given GKS basis J^, and then 
the basis on E^ is again GKS. 

One has the following ([8]) 

Proposition 2.8. If{E,fi) has a UOB which is GKS, then the GKSl in- 
equality is true. 

Proof. The previous statement just means that if we define a GKS function 
F as a function which may be written as F = '^j^Fifi, where (fi) are the 
elements of and Vi, Fi > 0, then if F and G are GKS functions, one has 



The statement is straightforward, since cxp(_ff) is again GKS, being the sum 
of a series with non-negative coefficients, and so F exp(ff ) is itself GKS. Since 
any GKS function has a non-negative integral, the conclusion follows. □ 

The GKS2 inequality is much harder. It has only be obtained in some 
restricted settings, like products of abelian groups, and when the basis comes 
from G^ for G elementary groups like dihedral groups, and some for other few 
groups. Nevertheless, in any example, one has both the GKS and the HGP 
property. 

There is no example of a GKS basis where the GKS2 inequality is not 
satisfied. But we may restrict ourselves to a simpler setting. 

Here is one conjecture that we had been unable to prove, and which mo- 
tivated most of the material of this section: 



Then 






Conjecture. // the UOB T has the GKS and the HGP property, then the 
GKS 2 inequality is true. 
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3 The hypergroup property in the infinite setting 

In the first section, we described the hypergroup property in the context of 
finite setS. In what follows, we consider a general probability space {E, £, /it), 
together with a basis T — {1 = /o, /i, . . . , /«, • • •)• Very soon, wc shall 
restrict ourselves to the case of a topological space (in fact an interval in the 
basic examples of Sections 4.3 and 5), where the functions /j of the basis will 
be continuous bounded functions. But some general properties may be stated 
in a more general context. 

3.1 Markov sequences associated with a UOB 

Let {E, a, general probability space. In this subsection, we shall ask {E, £) 
to be at least a "nice" measurable space in the context of measure theory. For 
us, it shall be enough to suppose that E is & separable complete metric space 
(a polish space) and that £ is the cr-algebra of its u-field. Then is 

separable. 

We suppose that some LP'{iJi) orthonormal basis T = (/o, . . . ,fn,- ■ ■) is 
given, with /o = 1. In what follows, we shall assume that this is a basis of the 
real Hilbert space although we may as well assume that the functions 

/„ may have complex values and be a basis of the complex Hilbert space. 
Such a unitary basis will be called a Unitary Orthonormal Basis (UOB) 
associated with the measure /U. 

We are interested in bounded linear operators K on L^(/x), for which the 
functions /„ are eigenvectors. They are uniquely determined by 



for some bounded sequence (A„). The central question we address here is to 
determine for which sequences (A„) one has 



for some Markov kernel k{x,dy) of probability measures on E. 

As before, we shall call such a sequence (A„) a Markov Sequence (MS in 
short) associated with the UOB J^. We shall say that the kernel k (or rather 
K with a slight abuse of notation) is associated with the MS (A„). 

In general, this is not an easy question, but as before the hypergroup 
property of the basis will be a way of describing all Markov sequences. 

Let us start with some basic remarks. 

First, for any Markov operator, K{fo) = /o, since /o is the constant func- 
tion, and therefore Aq = 1. 

Also, any such Markov operator is symmetric in i^(^), since we already 
know its spectral decomposition which is discrete and given by the basis J-'. 
That means that for any pair (/, g) of functions in (/z) , one has 



K{,fn) = Xnf, 
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KU){x)g{x) ^{dx) = / f{x)K{g){x) fi{dx). 



Therefore, the measure K (x , dy) fi{dx) is symmetric in {x,y). 

But any Markov operator is a contraction in L°°{fi), and any symmetric 
Markov kernel is a contraction in L^{n), since, for any / e L'^in), 

j \K{f)\ d,jL< j K{\f\) df, = j \f\K{fo) di^ = j I/I di^. 

Therefore, by interpolation, is a contraction in for any p > 1, and 

in particular in . 

We deduce from that that any MS (A„) satisfies 

Vn, |A„| < 1. 

Also, if (A„) and (pn) are MS's, with associated kernels K and Ki, for any 
e G [0,1], {eXn + {l-9)nn) is a MS, associated with the kernel 6'i^+(l-6')i^i. 

Therefore, the set of Markov Sequences is convex, and compact (for the 
product topology on M.^). This shows that describing all Markov sequences 
amounts to describe the extremal points of this convex set. 

Notice also that the set of all Markov sequences is stable under pointwise 
multiplication, which corresponds to the composition of operators. In other 
words, if is the set of extremal points in the compact set of Markov se- 
quences, and if {Xi{0)) is the MS associated with the point 6 G 0, then one 
has 

Xi{e)Xi{e') = [ Xi{ei) R{e,e',dei), 

Je 

for some probability kernel 6\ ddi) on the space 0. 

To determine that an operator if is a Markov operator starting from its 
spectral decomposition, we shall need the following proposition. 

Proposition 3.1. A bounded symmetric operator K on L'^in) is a Markov 
operator if and only if 

K{fo) = /o, / > ^ K{f) > 0. 

Proof. This is where we need the fact that the measure space {E,£,ii) is a 
nice space. The conditions on K are obviously necessary. To see the reverse, 
we apply the bi-measure theorem ([21], Page 129). We consider the map 

£x£^[0,l] : {A,B)^ J lAK{tB) dn{x). 

For any fixed B, this is a measure in A, and by symmetry, it is also a measure in 
B. Since our spaces are polish spaces, theses measures are tight, and therefore 
we may extend this operation into a measure fiK{dx, d,y) on the cr-algebra 
£ X £. The measure is symmetric, and any of its marginal is /i. 
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Then, we apply the measure decomposition theorem to write 

IJ.K{dx, dy) = K{x, dy)fi{dy). 
The kernel K{x, dy) is exactly the kernel we are looking for. □ 

3.2 The hypergroup property 

The GKS property is relatively easy to state in a general context, as soon as 
the functions fi of the basis are in L^{n), since then we may just ask that 

Vi, j, A;, / fififk dji > 0. 
Je 

But for the dual hypergroup property, one has to be a bit more cautious. 
In general, functions in are defined up to a set of //-measure 0. Therefore, 
the meaning of fi{x)/fi{xQ) is not so clear. In order to avoid difficulties, and 
since this shall correspond to the examples we are going to describe below, wc 
restrict ourselves to the following setting: is a compact separable Hausdorff 
space, and the functions fi are continuous on E. 

We may then set the following definition. 

Definition 3.1. We shall say that the UOB T has the hypergroup property 
(HGP in short) at some point xq G E, if, for any x G E, the operator defined 
on T by 

K (f) = M^f. 
fi(xo) 

is a Markov operator. 

In other words, we require the sequences [fi{x)/fi{xo)). to be Markov 
sequences. 

Observe that this implies that \fi\ is maximal at xo, since the eigenval- 
ues of a Markov operator must be bounded by 1. Since the functions /j are 
normalized in then for any z S N, one has |/i(xo)| > 1. 

As before, this definition is equivalent to the following 

Proposition 3.2. The UOB T has the HGP property at the point xq if and 
only if there exists a probability kernel K{x, y, dz) such that, for any i e N 



fi{x)fi{y) 
fi{xo) 



J 



= / fi{z) K{x,y,dz). 



Proof. We shall mainly use this in the obvious way: if there is a probability 
kernel K{x, y, dz) satisfying the hypothesis of the proposition, then the HGP 
property holds. In fact, if such a probability kernel K{x, y, dz) exists, for any 
X £ E, the Markov kernel K{x, y, dz) defines a Markov operator with Markov 
sequence {fi{x)/fi{xo)).. 
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For the reverse, if the HGP property holds, there exists for any x a Markov 
kernel kx{y,dz) which satisfies 



Si{x)fi{y) 



j fi{z) kx{y,dz). 



It remains to turn this family of kernels into a two parameters kernel 
K{x,y,dz). □ 

Sometimes, it is easier to see the HGP property in the reverse way. 

Proposition 3.3. If there exists a non-negative kernel k{x, dy, dz) such that 
for any i,j 

fi{x) 



I 



fiiy)fjiz) k{x,dy,dz) = Si 



''fiixo) 

then the HGP property holds at the point xq ■ 

Proof. Remark first that from our hypotheses, the kernel k{x,dy,dz) is a 
probability kernel (taking i = j = in the definition) . 

Now, the two marginals of the kernel k{x, dy, dz) are equal to /x, since 



/ 



fiiv) k{x,dy,dz) = 6oi 



which shows that those marginals and fi give the same integral to any /j, and 
therefore to any function. 

We may then decompose the kernel k{x, dy, dz) = ki{x, y, dz)fi{dy). 

Since the functions fi are bounded, we may consider the bounded functions 



Hi{x,y) = j fi{z) ki{x,y,dz). 



From the definition of fci and the hypothesis on k, it is straightforward to 
check that 

Hi{x,y)fk{x)fi{y) iJ,{dx)n(dy) = SuSik 



fi{xo)' 
and hence 

H.ix,y) = l^^0M. 

fi{xo) 

Therefore fci satisfies the hypotheses of Proposition 3.2 and the proof is 
completed. □ 

The representation of Markov sequences is then the same than in the 
previous section. 
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Theorem 3.1. // the UOB T has the hypergroup property, then any Markov 
sequence has the representation 



for some probability measure v on E. Moreover, the Markov sequences 



are the extremal Markov sequences. 

Proof. The proof is exactly similar to the finite case (see Theorem 2.1). □ 
Remark that the series 



does not converge in general. We shall see in the examples developed in the 
next section that the formal measure K{x, y, z)iJb{dz) is not absolutely contin- 
uous with respect to ji, and may have Dirac masses at some points. 

But we may still define a convolution structure from fi * fj = Sij fi / fi(xo) , 
which maps probability measures onto probability measures, and all Markov 
kernels associated with !F would be represented as K{f) = f * for some 
probability measure on E. We give no details here since this will not be used 
in the sequel. 

4 Sturm-Liouville bases and Achour-Trimeche's 
theorem 

4.1 The natural UOB associated with a measure on a compact 
interval 

In this section, we shall consider some natural infinite UOB coming from the 
spectral decomposition of Sturm-Liouville operators on a compact interval of 
the real line. 

Let us first describe the context. Consider a probability measure n{dx) = 
p{x)dx on some compact interval [a, b] C R. In what follows, wc shall assume 
for simplicity that p is smooth, bounded above and away from on [a, b]. The 
density p is associated with a canonical differential operator 






L{f){x) = f"{x) + ^{x)nx), 



The Hypergroup Property and Representation of Markov Kernels 



25 



which is symmetric in L'^{lj). We shall consider here L acting on functions 
on [a, b] with derivative at the boundaries a and h (Neumann boundary 
conditions). 

In this context, L is essentially self adjoint on the space of smooth functions 
with /'(a) = f'{b) = and there is an orthonormal basis 



of which is given by eigenvectors of L satisfying the boundary condi- 

tions. This means that there is an increasing sequence of real numbers 



From the standard theory of Sturm-Liouville operators, the eigenvalues Xi 
are non-negative and simple. Therefore, there is for any Aj a unique solution 
fi of the previous equation, which has norm 1 in and which satisfies 

/(a) > 0. We refer to any standard text book for details (see [15] or [40] for 
example) . 

This basis shall be called the canonical UOB associated with /i on [a, b]. 

The fact that we chose to deal with the Numen boundary conditions and 
not with the Dirichlct boundary conditions (/(a) = f{b) = 0) comes from the 
fact that we require the function 1 to be an eigenvector of the operator. 

It will be much more convenient in what follows, essentially for notations, 
to extend our functions by symmetry in a and b, (and the same for /i). In this 
way we may consider that we are working on functions on the real line, which 
are symmetric under x i-^ 2a — x, and are 2(6 — a)-periodic. 

The eigenvectors are perhaps not smooth then at the boundaries a and h, 
but they are at least (since they are solutions of the equation Lfi = —Xifi 
at the boundaries). 

The hypergroup property is stated at some point in [a, b]. In the finite case, 
we know at which point we may expect the hypergroup property to hold: this 
is a point of minimal mass. In the general case, such a reasoning does not 
hold, since one may choose a point with minimal density to a given reference 
measure, but this depends on this choice. 

Here, the basis is the sequence of eigenvectors of an elliptic second order 
differential operator L symmetric in i^(/x). In this setting, there is a natural 
distance associated with the operator L (in this precise example of Sturm- 
Liouville operators, this is the natural distance on R). In any example we 
know, the point Xq is minimal in the following sense 



:r = (i = /o,/i,... 



/nj • • •) 



= Ao < Ai < • • • < A„ < • • • 



such that 



A./„ ma) ^ .nib) = 0. 





We did not try to prove this in a more general context. However, it is not 
clear how the properties of the operator L must be reflected in the properties 
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of T to insure for example that the maximal values of the eigenvectors are 
attained at the same point, and that this point is of minimal mass in the sense 
of (3). 

4.2 Wave equations 

In this context, one has some other interpretation of the hypergroup property. 
On £> = [o, 6]^, we shall consider the following differential equation 

L,F{x,v) = LyF{x,y), (4) 

for a function F which has Neumann boundary conditions on the boundary 
of D. We shall say that such a function is a solution of the (modified) wave 
equation. 

We have to be careful here with the regularity assumption on the function 
F{x, y) that we require. We shall see later that given any smooth function f{x) 
at the level X = xo, with Neumann boundary conditions, there is exactly one 
smooth function F(x,y) on D which is solution of Equation (4) and satisfies 
F{x,yo) = f{x). 

In fact, if f{x) = aifi{x) is the orthogonal decomposition of /, then 

is a formal solution of the wave equation, since 

L^F ^Y,^^J^^M^)My) = LyF 

But we do not even know (for the moment) that this solution is such that 
LxF is in L^(/U ® ji). 

Therefore, we shall say that F is a weak L'^ solution of (4) if for any 
smooth function G{x, y) with Neumann boundary conditions on dD, one has 

j [(L, - Ly)G{x,y)]F{x,y) fi{dx)ii{dy) = 0. 

Since 

J L^{F)G iJL{dx)ii{dy) = J L^{G)F ii{dx)ii{dy) 

for any pair of smooth functions F and G satisfying the Neumann boundary 

conditions, then any ordinary solution is a weak one. 

Now, given any L'^{n) function f{x), the above construction produces a 
weak L'^{ii®jj) solution F(x, y) satisfying the wave equation (4), and we claim 
immediately that this solution is unique. In fact, writing the function F{x, y) 
as 
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F{x,y) = ^a,j/,(a;)/j(a;), 



and using the fact that the cigcnvahies arc simple, one may check that if F 
satisfies weakly (4) , then a^j =0 if i ^ j , from which we deduce our claim. 

Observe moreover that aoo = / F{x,y) ^{dx)^{dy), and that for almost 
every y^, F{x,yo) G L^{lj) and that 



It is not clear however that if F(x,yo) is smooth, then F{x,y) is smooth. 
This shall be done later at least when yo = a or yo = b. 

The link between solutions of the wave equation and Markov kernels is the 
following. 

If a Markov kernel is Hilbert-Schmidt (that is if its eigenvalues A, satisfy 
J2i < oo)) then it may be represented as a 



Therefore, there is a one-to-one correspondence between Hilbert-Schmidt 
Markov kernels and non-negative weak solutions of the wave equation 
which satisfy 



We then have the following 

Theorem 4.1. Assume that for any function f{x) on the interval [a, b] with 
Numen boundary conditions, there exists a unique solution Hf{x,y) of 
the wave equation (4) on [a, 6]^ such that H{x,ya) = f{x). Then, the HGP 
property holds at the point yo for the natural UOB associated with /i if and 
only if whenever / > one has Hf > on P. 

In other words, the hypergroup property is equivalent to the fact that the 
wave equation is positivity preserving. 

Proof. Assume first that the hypergroup property holds at the point yo. Take 

any smooth solution F{x,y) of the wave equation with F{x,yo) = f{x) > 0. 
Then, from what we just saw, one has 





where 



Hx,y) = ^\ifi{x)fi{y). 




F{x,y) = Ky{f){x) 



where Ky is the Markov kernel with eigenvalues fi{y)/fi{yo)- Therefore, 
F[x,y) is everywhere non-negative. 
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On the other hand, assume that any smooth solution of the wave equation 
which is non- negative on {y = yo} is non- negative everywhere. Consider the 
heat kernel 

Pt{x,z) =^exp{-\,t)fi{x)fi{z). 

We know that it is a smooth function on D, which is everywhere positive. 
Then, 

is the unique Lp' solution of the wave equation with Fc^^ix, yo) = PiX^, z). 

Therefore, this function is non-negative, and this shows that for any t > 0, 
the sequence 

exp(-A,t)A(£) 

is a Markov sequence. It remains to let i go to to get the result, since a limit 
of Markov sequences is a Markov sequence. □ 



4.3 Achour— Trimeche's theorem and wave equations 

In what follows, we consider the case of a symmetric interval [—6, b]. Then we 
have 

Theorem 4.2 (Achour— Trimeche). Let p he a log-concave and symmetric 
density on [—h,b]. Then, the natural UOB associated with fj, has the HGP 

property at the point — &. In this case, we may as well choose Xq = b. The 
same is true on any interval with any log-concave increasing density p. 

This result is one of the very few cases when one may produce hypergroup 
bases without any kind of group structure on the space E. We shall see in the 
next chapter that this property holds for Jacobi polynomials, but in this case, 
there are at least for the integer values of the coefficients some interpretations 
of the convolution which reflects the group action of some orthogonal group. 
There is absolutely no such interpretation in this context. 

In general, Achour-Trimeche's result is stated with a density p which 
vanishes on the boundary. Under the conditions usually stated in Achour- 
Trimeche's theorem, there arc then no difference between Neumann and 
Dirichlet boundary conditions. The series 



E 



h{x)h{y)h{z) , , . 



is absolutely continuous with respect to the measure p,, which is not the case 
here. 

We chose to present this result in the case where the density p is bounded 
from below because it seemed to us to be more natural. 
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Apparently, the proof of Achour-Trimeche's theorem had never been pub- 
lished. Wc found a mention of it in the reference book [13] and the result is 
announced in [1], with no proof. Most of the ideas presented here come from 
Achour's thesis. The idea follows a previous result of Chebli [16], which works 
on [0, oo) and is somehow simpler (It corresponds to the case of a concave 
decreasing density). 

Proof. To prove this result, we shall make use of the characterization of the 
hypergroup property in terms of the wave equation given in Theorem 4.1. We 
shall sec in the next paragraph that any smooth bounded function satisfying 
Neumann boundary conditions on [a, b] has a unique extension as a smooth 
solution of the wave equation (4), when xq is one of the boundary points 
(sec Paragraph 4.5). (This has nothing to do with the log-concavity of the 
measure or with the symmetry: this is just a consequence of the fact that 
\og{p) is smooth and bounded.) 

We first treat the case where the density p is log-concave symmetric. 

First we make use of the symmetry assumption. Then, any eigenvector of 
the operator L on [—6, 6] with Neumann boundary conditions is either even 
or odd, since fi{—x) is also an eigenvector with the same eigenvalue. 

Then, any solution F{x, y) of the wave equation (4), written as 

Fix,y) = ^aifi(x)fi{y) 

i 

is symmetric under the change (x,y) ^ {y,x) and under {x,y) (— .t, — y). 

Wc want to show that if F{x, —b) > 0, then F{x,y) > everywhere. For 
this, it is enough to show this on the domain Di = {x + y < 0, x > y}. 

Also, we may change F into F + epsilon for any epsilon > 0, and we 
arc thus reduced to prove that the result is true when the function / on the 
boundary is bounded below by some positive constant. 

Then, a point M e £>i, let Am be the triangle delimited by the lines 
X + y = c and x — y = c' passing through M and the line y = —c. Let M_ 
and M+ be the points of this triangle which lie on the line {y = —b}, M_ 
being the left point and M+ the right one (see Figure 1). Let F be a smooth 
solution of the wave equation (4), and let G{x,y) = F{x,y)p{x)p{y). 

2G{M) = G{M_) + G{M+) + ( G{s)a+{s) ds+ [ G{s)a_{s) ds 

JlM-M] JlMM+] 

(5) 

where 

and the integral J^j^j ^^^j H{s) ds and J^mm+] -^(^) '^^ denote the one dimen- 
sional integrals along the segments [M_M] and [MM+j against the (euclidean) 
length measure on those lines. 
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Fig. 1. Triangle Am- 



This formula relies on an integration by parts formula. Even though we 
shall use it only on domains like Am, it is perhaps of some interest to state it 
in general. So we set it as a lemma. We shall not give too much details here, 
since a more general formula will be derived in the next paragraph. 

Lemma 4.1. Let H a smooth function on D and Q d D a domain with a 

piece-wise boundary dfi. Then 

/ (Lx - Ly)H{x,y) p{x)p{y) dxdy = / WH n p{x)p{y) ds, 
Jn Jen 

where n = {nx,ny) denotes the exterior normal derivative of the domain and 

VH n = dxHux — dyHuy, 
ds designing the length measure on the boundary dO. 

We shall not prove this lemma. It is the analogue of the classical Stokes 
fornnila, where the elliptic operator A is replaced by the hyperbolic operator 
Lx — Ly, its invariant measure being p{x) p{y)dxdy . 

One may see this as a particular case of the general integration by parts 
formula 

j Hi{Lx - Ly)Hi p{x)p{y) dxdy = - j VHi VH2 p{x)p{y) dxdy, 

applied with Hi = 

Prom the previous formula, applied on J? = Am for a function F which is 
solution of the wave equation (4) and has normal derivative vanishing on the 
boundary [M_ , M_|_] , one has 
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- / + dyF)p{x)p{y) ds+ I id,F - dyF)p{x)p{ii) ds = 0. 

J[M-M] J[MM+] 

We may then perform a next integration by parts on both integrals to find 
G(M_) - G(M) + / G{s)a+{s) ds 

J[M_ M] 

+ G{M+) - G{M) + / G{s)a-{s) ds = 0, 

J[MM+] 

which gives (5). 

Under our assumptions of p, both a+ and a_ are non-negative on the 
subdomain Di: under the log-concavity assumption a_ is non- negative on 
{y ^ x}, and a{x) + a{y) = a{y) — a(— x) > if x -|- y < 0. 

Now, consider the smallest y & {—b, 0) such that there exists some point 
in {x,y) e Di with G{x,y) =0. On this point, we have 

2G(M) = > G(M_) + G{M+), 

which gives a contradiction. 

For the case where the density is log-concave increasing, we may use the 
same argument on the domain {x> y}, since we still have the solution of the 
wave equation symmetric under the change {x,y) ^ {y,x). Then we extend 
F by symmetry around the axes x = —b and x = b, and then by periodicity, 
into a function defined on R x [—6,6]. The same argument of integration by 
parts remains valid, and, by means of the symmetrization, the domain of 
integration (Mm+M_) that should be used is replaced by the same triangle 
Am as before, as shown in Figure 2. Then we use the fact that the function 
a is decreasing and non-negative. 
□ 

Notice that the second case (when a is non-negative) may be reduced to 
the first one if we extend p by symmetry around 6, into a log-concave function 
on the interval [—6,36], symmetric around the point x = b. (The function p 
may not be at the point x = b, but this causes no problem). Then one has 
to apply the previous result on symmetric functions on the interval [—6, 36]. 

4.4 Other representations of the solutions of the wave equation 

In this paragraph, we shall consider an operator L{f) = f" + a{x)f' in 
/ = [0,1], with Neumann boundary conditions, and look at different rep- 
resentations of the solutions of the wave equation (4) . 

We shall consider the probability measure p{x)dx in / which satisfies 
p'Ip = a- 

Let F(x. y) be a solution of the wave equation [L^ — Ly)F = on 
with Neumann boundary conditions. As before, it is easier to extend F to 
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Fig. 2. Triangle Am and (Mm+M_). 



by imposing symmetry conditions at the boundaries a; G Z or y e Z, and 
to extend a by imposing antisymmetry conditions on these lines (or if one 
prefers, symmetry conditions on p). 

We have seen before that such an equation has an integral representation 
(5). Our first task shall be to change it into a new one. 

As before, for M = {X.Y) with y > 0, we denote by Am the triangle 
delimited by the lines {x + y = X -{-Y),{y - x = Y - X) and {y = 0}. M+ 
and M_ denote the edges of this triangle which lie on the line {y = Q}. 

For S G Am, we denote by the unique point U on the interval 

[M_, M] such that S G [U, U+\ and the unique point U on [M, M+] such 
that S G [U-U] (see Figure 3). 

Recall that 

a+{x,y) = ^ {a{x) + a{y)), a-{x,y) = ^{a{y) - a{x)), 
and let R{x,y) = p{x)p{y). 

Proposition 4.1. If a continuous function G satisfies (5), then the function 

H{x,y) = ^tM= 
yR{x,y) 

satisfies 

2H{M) = H{M_) + H{M+) 

+ [ H{S)ao{M,S) dS+ [ H{S)a{M,S) dS, 

i[M_,M+] J Am 



The Hypergroup Property and Representation of Markov Kernels 33 

where 

a{M, S) = IVW) ( ^0^a.{S) + ^^^^^a+{S)] , 



and 



a,{M, S) = ' + 



\^/R{S^ ^/R{S+ 




+1 

Fig. 3. S, S+,S-, et S- 
Proof. We first start by considering the function 




and we notice that, for S G [M_, Af], 

eM{S) = 1 + ^ / 9M{U)a+{U) dU, 

^ J[S,M] 

and similarly that, for S G [M, M+], 

9m{S) = 1 + \ [ eM{U)a-{U) dU. 

^ J[M,S] 

Then starting from Equation (5), we replace the term 
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l[M-M] 

by 



/ G{S)a+{S) dS 

J\M- Ml 



/ {l-9M{S))G{S)a+{S) dS+ I 9M{S)G{S)a+{S) dS. 



'[M_M] J[M_M] 

In the last integral, replace G{S) by 



i \G{M-) + G{S+) + [ G{U)a+{U) dU + [ G{U)a-{U) du\ . 

^ \ J[M_S] HSS^] J 



2 

Then we have 



/ eM{S)a+{S)[ f G{U)a+{U) du\ dS = 

f G{S)a+{S) I f 9M{U)a+{U) du] dS, 

i[M_M] \^!76[S,M1 J 

and this last expression cancels with 

/ {l-eM{S))G{S)a+{S)dS. 

J[M-M] 

We do the same computation on the other side [M, M_|_] , and we collect 
the results. Observe that the term 

/ eMiS)a+{S)G{S-) dS 

gives rise to one part of the integral 

/ H{U)ao{M, U) dU, 

J[M-M+] 

[ [ G{U)eM{S)a+{S)a-{U) dSdU 



while the term 



'[M_M] J[SS+] 

produces one part of the integral 



H{U)a{M, U) dU. 



□ 
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We shall see in what follows that one may find many other integral repre- 
sentations of the wave equation. 
One is the following 

Proposition 4.2. Let F be a solution oj the wave equation (4), and let tp be a 
smooth positive function on P satisfying the Neumann boundary conditions. 
Then, if we set G{x,y) = F{x,y)ijj{x,y)p{x)p{y), we have 



2G{M) = G{M+) + G{M_) + [ G{s)K+{s) ds 

i[M_M] 



/ G{s)K_{s) ds - [ G(s) 

JlMM4.] J Am 



'[MM+] J Am 

K+{x, y) = V2{d, + dy) log (vVp(^)p(y)) 



ds, 



K-{x,y) = V2{dy - d^) log {'ij)^fl^(xjp(]j) 



where 
and 



Proof. The proof is the same as before, but we have to consider the equation 
satisfied by G instead of the equation satisfied by F. It is perhaps easier to 
make the computations under a change of variables 

u+v u—v 
X = — 7=-, y = 



V2' " x/2 ' 
in which case the operator — Ly becomes 

We extend our functions by symmetry to the set {y < 0}, and then 
our functions become for the new variables symmetric under the symmetry 
(u, v) I— > (w, u). Then the result is obtained through the integration on a square 
{uq < u < ui, uq < V < ui}. Since we shall not use these representations 
here, the details are left to the reader. □ 

As a consequence, if wc set U < V when U E Ay, and if there exists a 
positive function il) satisfying the Neumann boundary conditions with (Lx — 
Ly){ip) < and such that ^/p{x)p{y) t/j is increasing for this partial order, 
then any continuous solution of the wave equation (4) which is non-negative 
on {y = 0} is non-negative everywhere. In particular, if there is a solution of 
the wave equation which is increasing for this order, then the property holds. 

If wc arc looking for the hypergroup property at the point for the Neu- 
mann basis on [0, 1], a good candidate for the function ip in Proposition 4.2 
seems to be 
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tp{x,y) = 1 



where / is the (increasing) eigenvector associated with the first non-0 eigen- 
value, provided that |/(1)| < |/(0)|, which is a necessary condition for the 
hypergroup property to hold at 0. But we were unable to derive reasonable 
conditions on a which would insure that for this particular case the function 
^ p{x)p{y)ijj is increasing for the partial order on [0, 1]^. 

4.5 More about the solutions of the wave equation (4) 

As we saw in the previous section, there are many integral representations of 
the solutions of the wave equation on [0, 1]^. 
Most of them appear as 



where Vq (M, dS) is a continuous family of bounded measures whose support is 
the interval [M_ , M+] on the boundary {y = 0}, and Vi (M, dS) is a continuous 
family of bounded measures with support Am- 

In general, those representations lead to a unique representation 



for a continuous family of bounded measures with support in [M_, M+J. The 
crucial point is that in some situations the measure W{M, dS) may be positive 
even if Vi(M, dS) is not. In the case of Achour-Trimeche's theorem however, 
the measure W is positive only on some symmetric functions. 

To understand these representations, we shall consider a more general 
setting. 

Consider a separable compact Hausdorff space E, and two continuous fam- 
ilies Vo{M,ds) and Vi{M,ds) on E (two kernels). We shall identify such a 
family with the operator 



which maps the Banach space C{E) of continuous functions into itself. The 
identity operator corresponds to the kernel I{M, dS) = SuidS), and the com- 
position of kernels 



F{M) = / F{S) Vo{M,dS) + / F{S) Vi{M,dS) 




F ^ Vi{F){M) = [ F{y) Vi{M,dS) 



Je 




corresponds then to the operator composition. 
We set 
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/ \V{M,dS)\, 

J E 



- sup 

M Je 

which is the operator norm. 
Then we have 

Lemma 4.2. Consider some continuous function F e C{E) satisfies 

VM e E, F{M) = j F{S) Vo{M,dS) + j F{S) V^{M,dS). (6) 
// the series 

Ell^°"ll (7) 

n 

converges, then setting 

n 

one has 

F = W^(-)Vo{F). 

In particular, if Vq is supported by some closed subset Eq of E, and if the 
condition (7) is satisfied for Vi, then there exists a unique solution F to the 
equation (6) given the restriction Fq of F on Eq. 

Moreover, if the kernels Vq and V\ are non-negative and if Fq is non- 
negative, so is F. 

Proof. The proof is straightforward and is just the classical representation of 
(/ — Vi)"^ as X]„ V]^^"- Observe moreover that if V\ is non-negative, so is 
W\ = £iy\), and that the representation of the solution is then given by a 
non-negative kernel. □ 

In what follows, and to apply these lemmas, we shall consider the case 
where E \s a. compact subset of K x [0, oo[, and where, for some point M E, 
the measure Vb (M, dS) is supported by [M_ , M+] , Vi (M, dS) is supported 
by Am, and has a bounded density a{M, S) with respect to the Lebesgue 
measure on the product, in which case we write 



F{M)= [ F{S)Vo{M,dS)+ [ F{S)a{M,S) dS. 

i[M_,M+) Jam 



(8) 



Then wc have 



Proposition 4.3. Let k a uniform bound on \a{S,U)\, {S,U) G Am in (8). 
Then, on Am, we have 

||y0n|| <- K^I^mI" 



(n!)2 ' 

where \Am\ denotes the area of the triangle Am- 
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Proof. Recall the partial order {Si < S2) <s=^ e As^- 
Then one has 



yQn 



{M,dS)=lAM{S)an{M,S)dS, 



where 



an{M,S) = 



a{M,Sn-i)a{Sn-i 



Sn-2) ■ ■ ■ a{Si, S) dSi . . . dSn-i- 



S<Si<---<S, 



<M 



It is easy to see by induction that 



\[S,M]\ 



n 



an{M,S)\<K' 



.n 



where \ [S, M] \ denotes the area of the rectangle {U \ S <U <M). 
The conclusion follows easily from this estimate. □ 

Considering the representation of the solutions of the wave equation given 
in Proposition 4.1, we see that any continuous solution may be represented as 



where F(M, dS) has two Dirac masses at the points M_ and M+ and has 
a bounded density on (M_,M+). The smoothness of this density depends 
of course on the smoothness of the function a itself (and may be analyzed 
through the convergence of the series that we just described). For example, if 
a has k bounded derivative, then so has the density. 
We may observe the following. 

Corollary 4.1. Consider a solution F of Equation (6). Assume that Vq is 
non-negative and supported by Eq C E, and that V\ > V2, in the sense that 
V3 = Vi — V2 is a non-negative kernel. Suppose that V2 and V3 satisfy the 
growth condition (7) and moreover that £{y2) = W2 is non-negative. 

Then, if the restriction Fq of F on Eq is non-negative, then F is non- 
negative everywhere on E. 

Proof. Once again, this is straightforward. Setting W2 = (/ — ^2)"^, one has 



which is an equation of the same type, but with non-negative kernels. □ 




F{S) V{M, dS) 



F = W2 Vo{F) + W2Q V3{F) 



We may apply this for example for the solutions of (8) 
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Corollary 4.2. IfVo{M, dS) is non-negative and a{U, S) > —C, on S <U < 

M, where 

C= 



2|Z\m| 

and Ho is the first of the Bessel G function which is solution on (0, oo) of 
G" + — = -G, G(0) = 1, G'(0) = 0, 

X 

then any continuous solution of Equation (8) which is non-negative on the 

boundary {y = 0} is non-negative on Am- 

Proof. It is a simple application of the previous Corollary 4. 1 with V2 (M, dS) = 
-ClAM{S)dS. 

In this case, it is not hard to see that 

£{V2){M,dS) = 1aAS)F{\[S,M]\) dS, 

where | [S, M] \ denotes the Lebesgue measure of the rectangle [S, M] = 
{U I 5 < ?7 < M} and 



(-Ca;)" 



(n!)2 • 

The function F is the solution of 

xF" + F' = -CF, F(0) = 1, F'(0) = -C, 

which is related to the function G through the change of variable x = z^/AC. 
The function G is non-negative on [0,/xo) and this gives the result provided 
one observes that | [S, M] \ < 1^ . □ 

Remark. One may also observe that l//io is the fundamental eigenvalue 
of the Laplace operator on the unit ball of IR.^ with the Dirichlet boundary 
conditions, the function G(||a;||//io) being the corresponding eigenvector. 

All these considerations provide many criteria on the function a such that 
the associated Neumann basis on I has the hypergroup property at the left end 
point of the interval. In the next section, we shall deal with Gasper's theorem, 
where I = [0,n/2] and a{x) = atanx — /3cotx, with a > /3 > —1. The 
reader should check that no one of these criteria may apply on this example. 
Achour-Trimeche's theorem shows that the hypergroup property holds for 
this example in the symmetric case a = f3, even on any symmetric (around 
7r/4) subinterval of [0,7r/2]. But we do not even know for the moment if the 
hypergroup property holds in the general case on any symmetric subinterval 
of [0, 7r/2]. (Yet it is true for small subintervals and also provided that the 
parameters a and /3 belong to some specific domains that we shall not describe 
here) . 
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5 The case of Jacobi polynomials: Gasper's theorem 

Gasper's theorem states the hypergroup property for the family of Jacobi 
polynomials. The case of Jacobi polynomials may be considered as a special 
case of a Sturm-Liouvillc basis on [0,7r/2]. In this situation, both the GKS 
and the HGP property hold [24, 25, 26]. Actually, it is a unique situation for 
orthogonal polynomials, since they are the only ones, up to a linear change of 
variables, for which the HGP property holds (sec [18, 17, 19]) (under some mild 
extra condition on the support of the measure which represents the product 
formula). In the case of symmetric Jacobi polynomials (known as Gegenbauer 
or ultraspherical polynomials) , the HGP property may be seen as a particular 
example of Achour-Trimeche's theorem (although in this case the measure 
has a density which vanishes on the boundary). But in the general case, as 
wc already mentioned, none of the extension we gave of Achour Trimcchc's 
theorem covers this result. Even worse, we do not know if the HGP property 
holds for any symmetric subinterval of [0, 7r/2]. 

The Jacobi polynomials arc a quite universal object, since they are basi- 
cally the unique examples of a family of orthogonal polynomials which are 
also eigenvectors of Sturm-Liouville operators (together with their limiting 
cases the Hermite and Laguerrc polynomials, see [36]). On the other hand, 
for special values of the parameters, they may be considered as eigenvectors 
of rank-one symmetric compact spaces (here, with our notations, it is for the 
parameters (1,J5), (2,p), (4,p) and (p,p), with p G N). But for a wider range 
of parameters {p,q € N ) , they may be seen as eigenvectors of a Laplace oper- 
ator on a.p + q—1 dimensional sphere. The special case where p = q is much 
simpler, since then one may consider a p-dimcnsional sphere. 

In this section, after a short introduction on Jacobi polynomials and the 
statement of the hypergroup property for these polynomials, we present the 
simpler case of symmetric Jacobi polynomials, where the convolution struc- 
ture has a nice geometric interpretation for p G N. This interpretation is for 
example described in [12]. For the dissymmetric case, although the Jacobi 
polynomials still have a simple geometric interpretation too when the param- 
eters are integers, the convolution structure is far less obvious. 

5.1 Jacobi polynomials 

This polynomial family is defined for some positive parameters p and q as the 
family of orthogonal polynomials associated with the measure 

^lp,g{dx) =Cp,g{l-x)^{l + x)^dx 

on [—1,1], Cp^g being a normalizing constant such that Hp^g is a probability 
measure. 

These polynomials are also the eigenvectors of the operator 



The Hypergroup Property and Representation of Markov Kernels 



41 



on [—1, 1]. If Pj^'^ is the polynomial of degree fc, one has 

Remark. These polynomials are traditionally parametrized by a = and 
/3 = 2zl with a,/3 > -1, from [39] or [24, 25, 26]. 

If we change x = cos(2^), 9 G [0, |], then this operator is turned into 

LpM = \ inO) + {{q - 1) cot(^) - (p - 1) tan(^)) f{e)] . 

We see then that Jacobi polynomials is one example of a Neumann basis 
associated with a Sturm-Liouville operator (except that the density of the 
measure vanishes on the boundary points, for parameters larger than 2). We 
may also observe that the measure is log-concave as soon as p and q are in 

[l,oo). 

When p and q are integers, one may see the operator Lp^g as the action of 
some spherical laplacian on a quotient of the sphere. 

More explicitly, we set N = p + q. Let us denote by \X\ the euclidean 
norm of a point X in M^, and let be the unit sphere. We consider the 

Laplace operator Z\s on the unit sphere S'^~^ in M.^: this is the restriction to 
the sph(T(^ of the usual Laplace operator on acting on function which are 
defined in a neighborhood of the sphere and do not depend on the radius of 
the point. 

We parametrize S^~^ as 

where Xi e X2 G S'~^, and x G [—1, 1]. The action of on a function 
which depends only on x gives again a function of x, and we have 

As{h){x) = 4Lp,q{h){x). 

We shall say that such a function on the sphere which depends only on x 
(that is which depends only on the norm of the projection of X onto W) has 
the invariance SO{p) x SO{q), where the action of SO{p) x SO{q) is obtained 
by the action of the first component on Xi and of the second on X2. 

The measure /ip^g is the invariant measure for the operator Lp^q and the 
uniform measure is the invariant measure for the Laplace operator on the 
sphere. This shows that Hp^g is the image of the uniform measure on the 
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sphere (normalized as to be a probability measure) under the map X t-^ x of 
Formula (9). 

In fact, under this map X t-^ {Xi,X2,x), it is straightforward to sec that 
the uniform measure ap+q-i on §p+9-i is transformed into ap-i0aq-i '^ fip,q. 

Thanks to this remark, consider N > p and look at the projection tt{X) 
from g-^"^ onto the unit ball in MP. (That is the orthogonal projection when 
the sphere is imbedded into M^). If we set a; = 2 |7r(X)| — 1 and Xi — jT^xji ' 
we see that the image measure of crjv-i under X i— > {x, Xi) is Hp^N-p <S) cTp-i. 
This remark shall be used in Paragraph 5.4. 



5.2 Gasper's result 

Gasper proved the following product formula which gives the HGP property 
for Jacobi polynomials, applying Proposition 3.2. 

Theorem 5.1 (Gasper). Letp.q > and —1 < x,y < 1. Then 

• we have the following product form,ula: 

Vfc, ' = J P^iz) mp,q{x,y,dz) 

where mp^q{x,y,dz) is a Borel measure on [—1; 1]; 

• the measure rup q is positive (and then is a probability measure) if and only 

{P,q) ^{q>p}n{p>lorp + q> 4}; 

• moreover, if q > p > 1, mp^q{x,y,dz) is absolutely continuous with respect 
to Hpq, with density in L'^{iip q), so that 

with convergence of the sum for almost every z. 

The original Gasper's proof (see [25, 26]) consisted in the explicit computation 
of the sum Kp^q{x,y,z) using formulae on special functions like Bessel's and 
hypergeometric functions. There had been many other proofs of this property. 
For example, Koornwinder derived it in [35] from the addition formula of 
Jacobi polynomials and he found an other proof in [34], that we discuss next. 

Here we restrict ourself to prove the HGP property in the symmetric case 
{p = q) and in the case when q > p > 1. In the latter case, we follow a proof 
given by Koornwinder in [34]. However, his argument was based on an integral 
representation formula of the polynomials (our Lemma 5.2), whose proof, as 
we found in literature (see [3] together with [4]), relies on computational con- 
siderations on hypergeometric functions. In Section 5.4, we shall give a more 
geometric interpretation of this formula, at least when p and q are integers 
(it appears finally that the interpretation of Jacobi polynomials as harmonic 
functions was already known see [14, 23, 35] but it seems that it was not 
yet directly used to derive Koornwinder 's representation formula). 
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5.3 The special case of ultraspherical polynomials (p = g > 1) 

In the case when p = q, Jacobi polynomials are called ultraspherical polyno- 
mials. In this case, there is a much simpler representation of Lp^p when p gN* 
as the action of the Laplace operator on the sphere (and not on S^^~^ as 
before). 

When p is a positive integer, then the hypergroup property has a simple 

geometric interpretation, and thus the property is quite easy to establish. This 
easily extends to the case when p ^ N, by a simple extension of the formulae. 
This is what we are going to see in this paragraph. 

Consider a smooth function F : i-^- R which depends only on the first 
coordinate. To fix the ideas, let F{X) = f{X ■ ei), where ei is the first unit 
vector in M^'+^, and Y X denotes the standard scalar product in W~^^. Then, 
if Agp is the Laplace operator on S^, we have 

Asv = Lp,p{f)iX ■ ei). 

As before, the image measure of the uniform measure on the sphere 

through the map X i-^ x = X ■ ei is the invariant measure for Lp^p, that 
is /Xp^p. Moreover, we may parametrize S^'* = \ {ei, — ei} by 

X = (a;, Vl - x^Xi'^ , (10) 

where x & (—1, 1) is the first coordinate of the point x G EP C M^+^j and 
Xi e EP-'^. Through this map EP* ^ (-1, 1) x §p-\ the image measure of ap 
is /ip,p (g) CTp-i. 

From that, we see that if P^'^ is the A:-th ultraspherical polynomial, then 
P^'^{X.ei) is an eigenvector of A^p, with eigenvalue = —k{k +p — 1). 

Observe that for any point Y on the sphere, P^'^{Y.X) is again an eigen- 
vector on the sphere with the same eigenvalue. (This comes from the fact that 
the Laplace operator on the sphere is invariant under rotations.) 

Now, if we take two points Y and Z on the sphere, P|''^(y.X) and 
Pj^'^{Z.X) are two eigenvectors of A§p, with the same eigenvalue. Let us 
compute their scalar product in L'^{Ef) 

H{Y,Z)= f Pr{Y.X)Pr{Z.X)ap{dX). 

JSP 

Obviously, H(Y, Z) is a smooth function, taking values in [—1,1], and if R 
is any rotation, H{Y, Z) = H{RY, RZ). From this we see that H{Y, Z) = h{Y ■ 
Z), and we may write the function h in terms of ultraspherical polynomials 

h = ^arPr- 



We have 



44 Dominique Bakry and Nolwen Huet 

ar= [ h{x)Pr{x) Hp,p{dx) = I h{Y.Z)FP P{Y.Z) ap{dZ) 
Pr{Y.Z)Pl^^{Y.X)Pl'^{Z.X)ap{dX)aMZ). 



Using Fubini's theorem and the orthogonaUty of eigenvectors associated with 
different eigenvalues, we see that = unless r = k. We therefore see that 

/ P^'P{Y.X)Pr{Z.X) ap{dX) = a.PriY.Z). 

To compute Uk we choose Y = Z, from which we get, if we remember that 
the polynomials P^'^ have norm 1 in L , that 

akPrw = 1. 

Now, if we rewrite this formula for Y = ei and through the parametrization 
described above in (10), then we get, for Z = {z, y/1 — z'^Zi), 

/ pP'P{x)pP ''izx + Vl - zWl - x^X.Z,) Mp,p(dx)aj,_i(dX), 

while, for k ^ I, 

[ pP'P{x)Pl''P{zx + Vl-zW^-x^XiZi) Hp,p{dx)ap^i{dX) = 0. 
This may be rewritten as 

/ P^^^{x)Pf^^{zx + ^/l-z^^/l-xH) fip,p{dx)f,p.,,p.,idt) = Ski^ji^. 

Jx,t "k 

This last formula may be turned into an explicit representation 

PrWiy) kp{z,dx,dy) = 

for some probability kernel kp{z, dx, dy), which gives the hypergroup property 
thanks to Proposition 3.3. 

When p is not an integer, since we have an explicit representation of the 
kernel kp{z, dx, dy) = Kp{x, y, z)iJ.{dx)fi{dy), it is a simple verification to check 
that the function Kp{x,y,z) satisfies L^Kp = LyKp together with 

lim Kp{x,y,z) iipJdz) = 5x{dz), 

which is enough to get the HGP property at the point 1. 

Moreover, the convolution associated with this hypergroup structure is 
quite easy to understand when p is an integer. 
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Let us say that a probability measure /it on the sphere is zonal around 
X G S^* if it is invariant under any rotation R € SO{p+ 1) such that RX ~ X. 

Given any probability measure /x on [—1, 1], and any X G S^, we may lift 
H into a unique probability measure /t which is zonal around X such that the 
image measure of /t under the projection tt{Y) = Y ■ X from S^* onto [—1, 1] 
is fjb. 

Now, let us choose ei e and consider two probability measures v\ 
and V2 on [—1, 1]. We may lift vi into a probability measure vi on S^, which 
is zonal around ei. Then, we choose a random point in according to i/i. 
Then, given X, we consider the lift of U2 which is zonal around X and choose 
a random point Y according to this measure. Then, the resulting law of Y is 
zonal around ei, and we project this measure into a new measure vi * 1/2. It 
is an exercise to show that this convolution is the convolution associated with 
the hypergroup structure in this case. 

5.4 The case of dissymmetric Jacobi polynomials (q > p > 1) 

Although the dissymmetric Jacobi polynomials may be interpreted as eigen- 
vectors of the Laplace operator on the sphere 8*'+'^"^, it is far from trivial to 
prove the hypergroup property even in the case where p and q are integers. 
Nevertheless, the proof that we present below for completeness and which is 
due essentially to Koornwinder [34] has also some simple interpretation when 
p and q are integers in terms of harmonic analysis in R^+'J. 

Koornwinder's proof relies on two facts, given in the following Lemmas 5.1 
and 5.2. In what follows, and to lighten the notations, we remove the indices 
p and q from the definitions of the polynomials -P^"'- 

Lemma 5.1. (Bateman's formula) Let bk,r the the coefficients such that 




Then 



Pk{s)Pk{t) 



k 



{s + ty ( i + st \ 



r=0 



Lemma 5.2. (Koornwinder's formula) 



Pkjx) 



2{1 + x) - {1 - x){l + u) 
4 



[-1.1]=^ 



+ 



iV2Vl - + UV 

4 



Hp^q-p{du)fJbp-l^p-l{dv). 
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Before going further, let us show that this imphes the HGP property at 
the point xq — 1. In fact, wc shall use the characterization of the hypergroup 
property given by Proposition 3.2. 

For that, we replace in Bateman's formula of 5.1 the representation given 
by Koornwinder's formula. For this, we observe that, if {s,t) G (0, 1)^, then 



1 + 



> 1, 



s + t 

and therefore if we set x = {1 + st)/{s + t), we may replace by analytic con- 
tinuation — x"^ by va? — 1. 
Then, if we set 

, 2(1 + x)-{l- x){l -u) + V2Vx2 - iVTT^v 
i){s,t,u,v) = '- ^ ^ , 



one has 

Pk{s)Pk{t) 



From the definition of the coeflacients bk,r, we get then 
Pk{s)Pk{t) f Pk[{s + t)ip{s,t,u,v) - ij 



Pkiiy 



J Hp,g-p{du)np-i,p-i{dv). 



If we define mp^q{s,t,dz) to be the image measure of iip^q-p{du)iJp-i^p-i{dv) 
under the map 

{u, v) {s + t)tp{s, t, u, v) — 1, 

one gets 

Pk{s)Pk{t) /-p , , , ^ , , 
p^^^^ = J Pk{z) mp,q{s,t,dz), 

which is the announced result. 

Of course, one has to check that the image measure is indeed supported 
by [—1, 1], but this point is left to the reader. 

We now give the proof of Lemmas 5.1 and 5.2. As it shall turn out, they 
rely on elementary considerations on the interpretations of the operator Lp^g. 
For the moment, we restrict ourselves to the case where p and q are positive 
integers, and we shall interpret those formulae in term of the Laplace operator 
on MP+9. 

First observe that, given any function / on [—1, 1], we may lift this function 

on the sphere gP+'J"^ into a fimction which has the SO{p) x SO{q) invariance. 
Namely, using the parametrization of the sphere given in (9), we set 
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Notice that in this formula, 



where tti and ■K2 are the orthogonal projections on W and when the sphere 
is imbedded into Let us call U{f ) such a lift of a function from [—1,1] 

onto the sphere. 

Now, if Pfc is the Jacobi polynomial of degree k, the corresponding function 
Ffe = U{Pk) is an eigenvector of the Laplace operator on the sphere, and 
therefore the restriction to the sphere of a harmonic polynomial of degree 2k. 
Therefore, if we parametrize a point Z in M.p+'J by i? = \Z\^ and X = we 

may see that the function R^Fk(X) is harmonic in Rp+'. 

This may be seen in another way as we may write the Laplace operator in 
those coordinates 



where N = p + q and A§, is the Laplace operator on §^ ^. (It does not look 
as usual because of the change of r = |a;| into = R.) Since 



one may check directly that H{R, X) = R'^Fk{X) is a solution of AH = 0. 
In other words, the solutions of 



correspond to harmonic functions in W x which are radial in both com- 
ponents (bi-radial harmonic functions). If {X, Y) are the two component of a 
point in then this harmonic function is 



Proof. (Of Bateman's formula 5.1.) 

Let L = Lp^q. The function K{s, t) = Pk{s)Pk{t) is a solution of the wave 
equation 



In order to prove the assertion, which amounts to verify the identity of two 
polynomials, it is enough to check it on an open set. We shall choose to prove 
it on the set {s S (— 1, 1), t > 1}, on which the wave equation (Lg — Lt)K = 
becomes an elliptic equation. 

On the other hand, consider a solution G{R, x) of 







{Ls - Lt)K = 0. 
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■(f N d 1 ,^_„ 



and perform the change of variable 

1 + st 



R = s + t: X = 



t 



This equation becomes (L^ — Lt)G = 0. (We shall leave the computation to 
the reader, since it is just brute calculus.) 

This strange (and miraculous) change of variables may be much under- 
stood if we first operate a change of variables to reduce the leading terms to 
d1 + dy in both equations, and then observe that the transformation we made 
is conformal in M^, and thus preserves the leading terms. But we could find no 
simple geometric transformation, even in the case where p and q are positive 
integers, to understand this change of a bi-Jacobi equation into a bi-radial 
harmonic function. 

Therefore, the right-hand side of Bateman's formula is a solution of the 
wave equation [Lg — Lt)F = 0. The coefficients hk,r are computed in such a 
way that the two polynomials coincide out = 1. 

To see that they must coincide everywhere, it is enough to remark that if 
two polynomials A{s, t) and B{s, t) in (s, t) are solutions of the wave equation 
which coincide ont= 1, they coincide everywhere. Indeed, we may write 

k k 
A{s,t) = Y,arPr{s)Pr{t), B{s,t) = ^ 6^Pfe(s)Pfc(t), 
r=0 r=0 

and identifying the values in t = 1 produces Qr = W, r = 0, . . . , k. □ 

We now turn to the proof of Koornwinder's formula. 

We begin with a lemma. Here, we shall use for the first time that q> p. 

Lemma 5.3. Consider a bi-radial analytic function H on M^'+', that is to say 
H{X, Y) = h{ \X\^ , \Y\^ ) where h is an analytic function onM?. If moreover 
H is harmonic on 11$.*'+'^, then it holds 

H{X,Y)= [ H{X + m{RY),0) iy{dR) (11) 

JSO{q) 

where v{dR) is the Haar measure on the group SO{q), and w is the orthogonal 
projection from onto MP. 

As a consequence, if H{X,Y) = h{\Xf , \ Y\'^) and g{x) = h{x,Q), then 

h{\xf,\Yf) = 

J ^^^g(^\xf-^^\Y\^+iV2\X\\Y\VTT^s^ dnp,g-p{u)dnp-i,p^i{s). 

(12) 
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In practice, we shall just apply this lemma with polynomials functions h. 

Proof. (Of Lemma 5.3.) 

The proof comes from the following remark. Wc observe that if F is any 
analytic radial function in W^, namely F{X) = f[\X\^) where / is real ana- 
lytic, then F{X + iY) is a solution in Rp x MP of AxF + AyF = 0, that is 
to say that this function is harmonic in M^p. This is clear if we consider that 
F{X + y) is a solution of AxF = AyF. 

Remark here that in this analytic continuation, we consider some functions 
f{\X + iYf), where 

\X + lYf = \xf - \Yf + 2iX ■ Y. 

This is not the norm of X + lY considered as a point in C^. 

Then, F{X + m{Y)) is harmonic in Rp+'J, since the projection of the 
Laplace operator on W is the Laplace operator on MP. Hence, for any el- 
ement R e SO{q), F{X + 'nr{RY)) is harmonic in Mp+' since the Laplace 
operator on is invariant under rotations. 

Prom this, we see that 

H{X,Y)= I H{X + m{RY),0) i^idR) 

JsO{q) 

is harmonic. Observe also that it is bi-radial. It is obviously radial in Y , since 
we averaged using the Haar measure on SO{q). To see that it is radial in X, 
we just observe that, if R\ G SO{p), one has 

H{RiX + m{Z),0) =H{X + iR^\{Z), O) 

since H{X, 0) is radial. Moreover, for any Ri G SO{p), there exists R2 G 
SO{q) such that RnriZ) = tt{R2Z). 

Now, let us remark that we will get (12) from (11) just by expliciting 
the latter formula. For that, we write X = |X|ei, where e\ G S^*"^ and 

tt{RY) = \Y\ ^/^Yl, where Yi G Sp-\ 

Wc know that if R is chosen according to the Haar measure on SO{q), the 
law of -RjYj uniform on S'^ , and therefore, writing 7r(i?yp^) = \Y\ \J'^^-Yi, 
the law of (x, yi) is iJLq-p^p{dx) ® ap-i{dYi) , as we saw at the end of paragraph 
5.1. Therefore, if we set s = ei-Fi, the law of (u, s) is iJ.p^q-p{du)^iJ,p-i^p-i{ds). 

Then, 

\X + iTr{RY)f = \xf - \Yf + i\/2 \X\ \Y\ y/T+^ei ■ Yi , 

and 

/ g({\X + i7r{RY)f) v{dR) = 

JSO(q) ^ ' 

[ \Y\^+iV2\X\ \Y\ y/T+^s] iip,g-p{dx)fip-i,p-i{ds). 
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To finish the proof of the first formula (11), we observe that the two mem- 
bers of (11) coincide on Y = 0. On the other hand, the explicit formulation 
given in (12) shows that if g is analytic, then the right-hand side in (11) is also 
analytic in (|X|^ , |1^|^). Indeed, if we observe that the measure iip-i^p-i{ds) 
is symmetric, all odd powers of |F| |X| in the polynomial extension of 

1^1' - ^ \Yf + iV2 \X\ \Y\ VTTr 



will disappear through integration. We are therefore left with a series in 
{\X\\\Yf). 



It remains to see that two analytic harmonic bi-radial functions which 
coincide on F = coincide everywhere. Let , \yf) an analytic bi-radial 

harmonic function on The function ft, is a solution of 

(xa^ + |5, + 2/a2 + |a,)/i = o. 

Then, we see that if we write the expansion 
one has 

(n+l)(n+p/2) 



{m + l){m + q/2)' 

This shows that as soon as one knows (a„_o)) one knows h. This completes 
the proof of Lemma 5.3. □ 

Proof. (Of Koornwinder's formula (5.2), for p and q integers. ) 

In the case where p and q are non-negative integer, it turns out that it once 
again relies on properties of the harmonic functions in the Euclidean space. 

First, we lift both members on Mp+* and then we multiply them by R'^, 
where R = \X\^ + \Yf. 

As we have seen before, the function 

{\xf + \Yf)''Pk 




is a bi-radial harmonic function, which is a polynomial in \Xf and l^l^. 
It remains to apply Lemma 5.3 to conclude the proof. □ 

If we want to extend the proof of Koornwinder's formula when p and q are 
no longer integers, then we just have to observe that, setting S = \Xf and 

2 

T =\Y\ , we used the fact that the fact that 



S-T 
S + T 



The Hypergroup Property and Representation of Markov Kernels 51 
is a solution on (0, oo)^ of 

(^Sdl + ^ds + Td^ + I^t) H = 0, 
and that for any analytic function /, the function 

K{S,T) = J f - ^~^^ T + \V2VST^1 + us^ iip^q_p{dx)iip_i^p_i{ds) 

is also a solution of the same equation (but this time, one has to compute 

that by brute force!). 

Remarks. 1. The proof of Koornwindcr's formula gives a representation, 
for analytic functions, of solutions H{S, T) of LH = 0, where 

L = Sdl + ^ds+Td^T + \dT, 

in terms of the boimdary values H{S,0). This is some kind of Poisson 
formula. In such a formula, one has (at least for bounded functions) 

H{x,y)=E,,y{H{XT,YT)), 

where {Xg, Yg) is the diffusion with generator L, and T the hitting time 
of the boundary. 

Here, at least when q > 2, the boundary is polar and the set {y = 0} is 
never attained. But the representation is given here through a complex 
variable (and of course our functions are unbounded) . So it happens "as 
if" the process is willing to hit the boundary, provided one allows complex 
values (we do not know which is the meaning of that, of course). But it 
is certainly worth looking for more general integral representations of this 
type, with complex values on polar sets. 

2. Of course, when q converges to p, the measure iJ,p^q-p{ds) converges to the 
Dirac mass at the point 1, and Koornwinder's formula of Lemma 5.2 gives 
for the ultraspherical polynomials 



One may check directly this formula when p e N, through a much simpler 

argument, using harmonic functions in W^'^ instead of harmonic functions 
in M.'^P. This time, one has to extend the polynomial P^'^ into 

i^.(X) = |X|'=Pr(^-ei), 
where ei is any point of the unit sphere. 



.VT 



f^p—l,p—l 



{dv). 



(13) 
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